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ABSTRACT 



All known adaptive beam formers utilize some form of 
automatic minimization of the mean square error. High 
adaptation rates though, exhibit a signal cancellation 
phenomenon leading to self- jamming by the adaptive antennas. 
This effect results from adaptive interaction between signal 
and interference (i.e., jammer) inputs simultaneously 
received by an adaptive antenna. This research investigates 
various existing ways of adaptive beam forming for noise 
cancelling, and signal enhancement from simple Adaptive 
Noise Cancellers to Hard Constraint Adaptive Beamformers. 
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I. INTRODUCTION 



Adaptive antennas have been under development in variou 
forms during the past two decades, having thus far proven 
themselves capaole of rejecting various jamming signals. 
Most high performance radar and communication systems be: ' ; 
designed to work in jamming env 1 ro nmen t s currently 
incorporate various forms of adaptive antennas. Tie 
concurrent development of spread spectrum techniques •; ich 
the adaptive antennas provides a set of technologies for 2 
resistant systems. Both of the above technologies are 
compatible and frequently are used in the same system. To-: 
adaptive antenna is to attenuate the strong jamming signals 
as they appear at the front end of the receiver; then spree 
spectrum techniques are used to neutralize a large number 0 
weak signals that may not be eliminated by the adaptive 
antenna, and to recover the desired signals. 

The way an adaptive antenna works is to pass the signal 
contaminated with noise through a filter that tends to 
suppress the noise while leaving the signal unchanged. 

These kind of filters, known as adaptive filters, have the 
ability to adjust their own parameters automatically and 
their design requires little or no prior knowledge of the 
signal or noise character istics . 
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Adaptive noise cancelling is a variation of optimal 
filtering that uses an auxiliary, or reference input. This 
input is then subtracted from the primary input, which is 
composed of both the desired signal and the undesired noise. 
This procedure of filtering and subtraction is controlled by 
an adaptive process, giving noise reduction which, depending 
on the particular application, can result in significant 
noise cancellation. 

The design of an adaptive antenna has to retain signal 
preservation. Unfortunately some adaptive beam f or nets 
use do not perform well in certain environments wnere 
correlated -signal conditions lead to partial or total 
cancellation of the desired signal witnin the antenna. 

In that way, although our primary intention is tne 
cancellation of all the undesired interferences, we create a 
self-gamming of the adaptive arrays by the so-called 
phenomenon of signal cancellation resulting from the 
adaptive interaction between the signal and the jammer which 
simultaneously are received by the adaptive antenna's 
elements . 
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ADAPTIVE NOISE CANCELLER 



II . 



A. SASIC CONCEPT 

The oasic block configuration of an Adaptive Noise 
Canceller is shown in Figure 2.1. A signal, s, is 



_i a.: o.il . w 



ted over a channel to a sensor wnicn also rec-r 



noise, n , uncorrelated witn tne signal. 



c omo i 



second sensor receives a noise, n n , also uncorrelated w 
tne signal, out correlated in some unknown way with tne 
noise, n . This second sensor provides the so-called 



referenc 


e input of tne 


cancel leu . 


Noise, n , , is 


prod uce 


the output, y, 


tnat is as 


close as possio 


replica 


of n . Out out 
o 


, y, is now 


subtracted from 



primary input to give the final system's output. 



z = s + n - y 



(2.1) 



The reference input, n^, is processed by an adaptive 
filter which automatically adjusts its own impulse response 
This adjustment is accomplished through an algorithm that 
responds to an error signal dependent on the filter's 
output. In this way, using the proper adaptive algorithm, 
the filter can operate under changing conditions which 
minimize the error signal. 
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Figure 2 . 1 . Basic Configuration of the 



The error signal, e , is feeding the system ' s output oac> 
to t-he adaptive filter in order to adjust the filter through 
a LMS (Least Mean Squares) adaptive algorithm and minimizes 
the total system output power. In other words, in an 
adaptive noise canceller the output forms an error signal 
for the adaptive process. 

B. ANALYTICAL MODEL— NOTCH FILTER RESPONSE 

In order to analyze the adaptive filter, we coo os 
pure sinusoidal for reference input 



Glover [Ref. 1] has shown that under this condition etc- 
dashed box of Figure 2.2(a) can be approximated by Figure 
2.2(b) . 

The conventional way of eliminating such sinusoidal 
interferences is through the use of a notch filter. Figure 
2.3 shows a single frequency noise canceller, with two 
adaptive weights, which are updated through the equations: 



x . = C cos ( j t + $ ) 



( 2.2 




(2.3) 




The sampled reference inputs are: 
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Figure 2.2. 



Analytical Model of the A.tT.C. 
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SAMPLING FRLQ. Q = ~ RAD/SLC 



(2.5) 



= C COS(u q ]T + ?) 

x 2J = C sin(ii Q ]T + 5) (2.5) 

Widrow (Ref. 2] and [Ref. 3] has snown that the single 

frequency noise canceller has the properties of a notch 

filter 1 at the reference frequency, j , having the tran^i-'. 

o 

function : 



Y (z) /X (z) = H (z) 



whe r e 



H (z) 




z 2 - 2z cos (2 n mi 8 ^ ) + 1 

- jC 2 )z cos (2’rjh ^ ) + 1 - 2 

o 



The zeros of the above transfer function are located - 
the z-plane at 



_ ±J2 mu a 

Z = e O 



( 2 . 8 ) 



and are inside the unit circle. 

The poles, also inside the unit circle, are located at 



z = (1 - uc 2 ) cos ( 2 iru a ^ ) ± 

j [ (1 - 2 uc 2 ) - (1 - uc 2 ) cos 2 ( 2 ttco o ) ] 1/2 (2.9) 

^Detailed explanation is given in Appendix A. 
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Figure 2.4 indicates the location of poles and zeros m 
the z-plane. Since the zeros lie inside the unit circle, 
the depth of the notch in the transfer function is infinite 
at the interference frequency, w . The sharpness of tne 
notch is determined by the closeness of the poles to the 
zeros . 

Figures 2.5, 2.6 and 2.7 are magnitude plots of z 

transfer function H ( z j using as C tne value of one, 
different /alues of tne adaptation cons tan- -. 

C. NOISE CANCELLATION' IN THE NARR0W3 AND CASE 

Consider tne case in wnich both signal ana jammer 
interference) are pure sinusoids of known frequencies, f ,, , 

f j , where 



Figure 2.3 shows the frequency spectra of both noise an 
gammer, and Figure 2.9 indicates the same spectra after 
being filtered by a notch filter with cutoff frequency the 
same as the one of the jammer, f^. 

Comparing the two preceding Figures 2.3 and 2.9, we see 
how effective the notch filter is in eliminating sinusoidal 
interferences. Of course, pure sinusoidal interferences ar 
seldom found in nature, and still the main disadvantage is 
that noise and jammer have to be of different frequencies. 
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Figure 2.5. Frequency Response ol n Not:rh !• i M or (u 0.001) 
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Fiqure 2.6. Frequency Response ol <i Notch Fill cm ( n 



WOVRPI = W/P 1 , MO. I i riUNi) 



C 



02 *7 




06 '0 



SL -0 



09 '0 



ShO 

ObH 



Q£ ’0 



St -0 



00 ‘(T 



21 



Fiqure 2.7. Frequency Response 



rfcxmmer 
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Ficjure 2.8. Frequency Response ol S i nuso i <l,i I S i < j na ) - Jnmme r 
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Figure 2.9. Frequency Response oi Two Fillon**! Sinusoids 



III. 



FROST BASED ADAPTIVE 3EAMF0RMER 



A. CONSTRAINED LEAST MEAN SQUARES ALGORITHM 

Frost [Ref. 4] proposed an algorithm, the L.M.S. 
algorithm, for adjusting an array of sensors in real tire r 
respond to a desired signal wmie discriminating a j - i .1 s -. 
noises. Figure 3.1 illustrates the basic set-up for 
implementing an L.M.S. algorithm. It is assumes on - o to- 
desired signal arrives in plane waves from a c man 
direction called the look direction. 

The algorithm iteratively adapts weights m ard-r 
minimize noise power at the array output while ma mt-. - 1 ’ n j 
chosen frequency response in the look direction. Require- 
ments for the algorithm is a priori knowledge of t.ne 
direction of arrival and the frequency band of interest. 
During the adaptive process, the algorithm progressively 
learns the statistics of noise arriving from all directions 
except the look direction. 

A major advantage of the constrained L.M.S. algorithm i 
that it has a self-correcting feature, permitting it to 
operate for arbitrarily long periods of time in a digital 
computer implementation without deviating from its 
constraints oecause of cumulative roundoff or truncation 
errors . 
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Figure 3.1. B£isic t'on t iqu r u l i on ol 



Noise arriving from the look direction may be filtered 
out by a -suitable choice of the frequency response 
characteristics in that direction or by some external means. 

B. BASIC PRINCIPLE OF THE CONSTRAINTS 

L.M.S. Algorithm maintains a chosen frequency response 

in the look direction while minimizing output noise power 

because of a simple relation oetween the look direction 

frequency response and the weights in the array of Fig ire 

3.1. Assuming that the look direction is chosen 

perpendicular to tne line of sensors, identical signal 

components, arriving on a plane wavefront parallel to tne 

line of sensors, appear at the first taps simultaneously aid 

o 

parade in parallel down the tapped delay lines" roll-twin; 
each sensor. However, noise waveforms arriving fron other 
than the look direction will not, generally, produce equal 
voltage components on any given vertical column of taps. 

The voltages, signal plus noise, at each tap are multiplied 
by the tap weights and added to form the array output. 

Thus, as far as the signal is concerned, the array processer 
is equivalent to a single tapped delay line in which each 
weight is equal to the sum of the weights in the 
corresponding vertical column of the processor, as indicated 
in Figure 3.2. These summation weights in the equivalent 

2 

Appendix B describes the tapped delay line filter. 
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tapped delay line ,nust be selected so as to give the desired 
frequency response characteristics in the look direction. 

If the look direction is chosen to be other chan rhe one 
perpendicular to the line of sensors, -hen the array can oe 
steered either mechanically, or electrically, by the 
addition of steering time delays placed immediately afcer 



eacn sensor . 

Mow consider the array proole.m of Figure 3.2. 
In vector notation the sample tap voltages ire 
id en 1 1 f 1 ed as 

X T (3) = [x 1 (g), x 2 (j), ..., x j ( ] ) ] 

and the sampled weights as 


: . . 


W T (j) = [w 1 (j), w 2 (j), ..., ^ kJ (g)] 


( 3 . 2 J 


where g represents the sample number. The filter ou‘ 
are summed to form the beamformer output 


:puts 


y ( j ) = W T (j) • X(j) = X T (g) • W(j) 


(3.3) 


If we denote by C the constraint matrix, defined as 





c [ c ^ , c 2 / • • • t c j ] 


(3.4) 
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a KJ xJ dimensional matrix, that conveys the beamform-er 
structural info-rmation needed to implement the constraints, 



then the constrained L.M.S. problem may be stated as 

Minimize E[y^(j)] = E [ W T • X ( j ) X T ( j ) ♦ W] (3.5) 



subject to 


C T W = .-/L (3.5/ 


wn ere is :ne 


looK-d i rect ion response vector defined as 



_ ~ [£-^r ^2' ( 5 . 

The weights are updated using the constrained L.M.S. 
algorithm, which for the case of Figure 3.2 (assuming two 
sensors) , can be written as 



w x ( : +i) 


= (l/2)[w 1 (j) - uy(j)x 1 (j)] 




- (1/2 ) [ w 2 ( j ) - t>y ( j ) x 2 ( j ) ] + (f x /2) (3.3) 



w 2 ( j +1 ) 


= ( 1/2 ) ( w 2 ( j ) - uy(j)x 2 (j)] 




- (l/2)[w x (j) - uy ( j ) x x ( j ) ] + (f^/2) (3.9) 
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w 3 ( j + 1 ) = ( 1 / 2 ) [ w 3 ( ] ) - uy(j)x 3 (j)] 

- ( 1 / 2 ) [ w 4 ( j ) - wy(])x 4 (])] + ( f 2 / 2 ) (3.10) 

w 4 (j+l) = ( 1 / 2 ) [ w 4 ( j ) - uy(j)x 4 (j) 

- ( 1 / 2 ) [ w 3 ( 3 ) - u y ( j ) x 3 ( j ) ] + (f/ 2 ) (3.11) 

The performance that a FROST oased adaptive beam for ner 
exhibits is unity gain and flat frequency response in the 
look direction. A deep null must oe placed at the jammer's 
frequency. Figure 3.3 is a plot of the signal direction 
frequency response, and Figure 3.4 is a plot of the jammer- 
direction frequency response. Figure 3.5 indicates the 
beampattern of the FROST beamformer, confirming that a deep 
null has been placed on the jammer. 

C. NOISE CANCELLATION IN THE WIDEBAND CASE 

The first case for -examination is the one of a 
sinusoidal jammer of frequency, f Jf and a broadband signal 
of center frequency, f , where 

o 

fj = f s (3.12) 
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Figure 3.3. S ignu 1 -T)i red i.oti Freguenry 
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Fiqure 3.4. Jammer -Hi red i on \'\ o* juenrv i«* soon so 
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Figure 3.6 represents the spectrum of the signal, since 
Figure 3.7 indicates the spectrum of the stronger sinusoidal 
gammer. Figures 3.3, 3.9 and 3.10 represent results of 

passing the above described signal and noise through notch 
filters, each one employing a different adaptation constant. 

Comparing the preceding three figures we see now 
important the adaptation constant, u, is and now it can lead 
to signal aliasing by choosing it to have bigger than proper 
mag n i tude . 
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Fiqure 3.6. F r eq uericy Kesuonsu 
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Pig u re 



IV. DUVALL BASED ADAPTIVE 3 SA.MFORMER 



A. BASIC CONFIGURATION 

We have already illustrated some relatively simple 
s igna 1/j atnmer scenarios and have seen than nor only the 
Adaptive Noise Canceller (A.N.C.) , but even the .tost ng.i 
constrained beam formers (Frost A.B.F.) can fail to preset /e 
the desired signal from distortion. 

In this chapter we will review and describe a n e n e - 
solution to the noise cancellation problem and see r.z .. we 
can preserve the desired signal at the price of some 
increase in beamformer complexity. 

The basic configuration of the DUVALL based been ; /c; 

[Ref. 5] is illustrated in Figure 4.1. 

Two observations were made by Duvall useful in 
developing the Composite Beamformer (C.B.F.) . The first 
observation was that interaction between the desired signal 
and the jammer is the root for the signal cancellation 
phenomenon. Frost experimentally and analytically has shown 
that the presence of both signal and jammer energy is a 
prerequisite for signal cancellation, and the output signal 
containing correlated signal and noise waveforms is the 
phenomenon that makes the signal distortion occur. 

The second observation was that the signal plays no role 
in the optimum (i.e., Wiener) solution calculation in a 
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*Nulls the desired siqnal 



perfectly steered Frost A.B.F.. In other words, the look- 
direction response is determined exclusively by the hard- 
constraints and not by the characteristics of the desired 
signal . 

These two above observations were addressed for 
exclusion of the signal from the beam former . 

Reconsidering uhe whole structure of the C.B.F., we see 
that the key elements are: 

1. An Augmented Array . Tne array has been augments’: so 
that subarrays consisting of multiple elements appear 
in place of the individual elements of one original, 
array . 

2. A P reprocessor . A preprocessor operates upon one 
received signals from the augmented array oo genet- oe 
an environment that is free of desired-signal content. 

3. An Adaptive 3eamformer . The adaptive beam, former is 
required that can be constrained to control one loot- 
direction response, while nulling all gamming signals. 

4 . A Slaved Non-Adaptive Beamformer . The ad a p t i v e 
beamformer derives weights that are copied to one 
slave beamformer, which has the same structure as one 
adaptive beamformer but is connected directly to 
selected antenna elements and is used to implement tne 
computed solution and recover the desired^ signal. 

Three of the above key elements, the array, the 
preprocessor, and the adaptive beamformer, afford consider- 
able flexibility. A variety of specific realizations are 
possible. The slaved beamformer design is inflexible in the 
sense that it mirrors the design of the adaptive beamformer. 



B. PRINCIPLES OF THE DUVALL C.B.F. 

The general C.B.F. structure of Figure 4.1 can be 
specialized to yield the C.B.F. shown in Figure 4.2. Here 
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Fiqu re 



the desired signal is assumed to be incident: from broadside 
so we can neglect beam steering. The preprocessor is 
realized by using two-element subarrays in a simple element; 
differencing scheme. The snaring of elements oecween 
subarrays provides very efficient use of elements; only a 
simple element is needed beyond chose ordinarily required 
for a comparable array. 

In order to prevent the desired signal from distortion 
(i.e., signal cancellation), the signal wnicn arri/es from, 
tne look-direction is excluded from the beam f o rmer .-/mere tie 
adaptive process takes place. Under the proposed senate, 
the adaptive process is only used to determine a set of 
weights. These weights are copied into a separate, 
identical processor which is used to form the output signal. 
The C.B.F. scheme can be in cooperation with any existing 
adaptive algorithm including the Frost one, as Figure 4.2 
add res ses . 

Due to the subtractive preprocessing, the look-direction 
signal does not appear at the Frost adaptive beamformer 
inputs, whose inputs are only composed of the jamming 
signal. The task of the Frost Algorithm is to place a null 
at the gammer direction through the determination of an 
appropriate set of weights and the signal has now no effect 
at all on the determination of those weights. By copying 
the weights the slaved processor places the main beam of the 
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antenna pattern to the constraints established by the look- 
direction and verifies that a null has already been placed 
in the jammer's direction. 

Figure 4.3 represents the Duvall C.B.F. in terms of 
phasor notation. The jammer components received by the 
antenna are indicated by a set of equal amplitude unifocm _y 
spaced phasors, Jq, J ^ , J-,, and The phasor inpits t 

the Frost beamformer are J, - J Q , - J^, J-j - J 9 , and 

j. - j,. They are also un i f orm-ampl i tude , equally-spaced , 

rx J 

and separated by the same angles as the received jammer 
components, Jq, J^, 3 2 > J 3 an d -^ 4 * 

Since the relative phase angles are the same in the 
slaved processor as in the Frost processor , the beam 
pattern notch is formed at the proper bearing angle. 

The phasor argument applies to a single jammer at a 
single frequency. Linearity and superposition apply and 
show that phase relations are preserved for multiple jammer 
and for broadband as well as narrowband signals. 

The uniform linear array provides an attractive 
structure for the C.B.F. because there is the option of 
element sharing between subarrays. A regular array 
structure is not, however, a prerequisite for the C.B.F.. 
The fundamental requirement is for phase matching between 
the master and the slave beam formers . Phase matching may b 
obtained for an arbitrary array geometry by augmenting each 







-3 



46 



Figure 4.3. DUVALL B . F . Humor Diagram 



original element to form identical subarrays at the e 
locations. Identical preprocessors may then oe used 
subarray responses with nulls in the selected look 
direction. The preprocessor outputs must be cophased 
the look direction and applied to the adaptive beainfo 
wniie the original element outputs are cophased and a 
to the slave. 
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V. 



SPATIAL DITHER ALGORITHMS 



A. GENERAL IDEA 

Spatial dither algorithms [Ref. 2] have been newly 
conceived for the purpose of applying locally controlled 
modulation to signals arriving at angles otner man me 
look-direction, while leaving inputs arriving from the !oc'- 
direction uncnanged . Tne technique that these algorimcs 
are employing is focused on a reduction of the mm-, p^w-i' 
density by spreading it spectrally. When me spatial dir 
is used with a conventional adaptive beam former , it red ices 
tne signal cancellation effects. 

Conceptually, a simple form of spatial dicner algor; mm 
is represented in Figure 5.1 and stated as the "3 4-mco 
plywood" approach. The elements of an antenna array may oe 
imagined to oe fixed to a piece of plywood wnich provides 
rigid insulating support, so that the entire array may be 
moved mechanically. In either one or two dimensions, the 
array is moved in directions which are orthogonal to the 
look-direction. Far-field ananations arriving from the 
look-direction will be undistorted by the mechanical motion, 
while emissions from off axis sources will be distorted by 
an unusal s’ni f t-o f- time-base form of modulation. 

The outputs of the antenna elements of Figure 5.1 could 
be applied to a time delay and sum ( nonadaptive) beamformer. 
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Figure 5.1. Spatial Dither Algorithm (3/4” Plywood Approach) 
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to a conventional adaptive beamformer, or to a Duvall 
adaptive beam former. Spatial dither could be beneficial 
each case. By reducing jammer power coherence, some an: 
jam protection is provided without adaptive beam fo r:n ing , 
additional anti-jam protection is provided with adaptive 
beamforming. Reduction of signal cancel lac ion efface 
Frost beam former can be obtained by using special ditner 
preprocessing. Breakup of jammer signal structure is a 
possiole form of signal preprocessing applicaole co all 
cypes of adaptive and nonadaptive beam fo rimer s . 
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VI . 



CONCLUSIONS 



Adaptive Noise Cancelling is a method of optimal 
filtering that can be applied whenever a suitable reference 
input is available. The principle advantages of the ri;id 
that an A.N.C. employs are its adaptive capaoility, its 1 ;u 
output noise, and its low signal distortion. The adap-: 
capability allows the processing of various signals ,;noso 
properties are unknown. Output noise and signal distort: u 
are generally lower than has been achieved with the v a r Lo u - 
conventional optimal filtering configurations. 

The FROST beam former has proven itself good _;i regent 
various interference signals and enhancing the signal of 
interest for several cases. Interaction, thougn, oetween 
the desired signal and the various interference signals car. 
lead to partial, or total, cancellation of the desired 
signal within the adaptive beamformer. 

The DUVALL composite beamformer gives a solution to the 
problem of signal cancellation, at the price of increased 
complexity in the beamformer structure and implementation. 
Under that scheme, in order to prevent signal cancellation, 
the useful signal is excluded from the beamformer in which 
the adaptive process now takes place to derive a set of 
weights. These weights are now copied into a separate, 
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identical "slaved" processor used to form the output signal. 
The Duvall approach to the problem is appealing but since it 
is a new development, possible limitations of its 
performance have not yet been assessed. 

A WIDROW solution to the problem involves spatially, or 
electronically, moving the receiving array to nodulate 
emanations received off the look-direction wirnour 
distorting useful signals in the look-direction. Tnis 
approach, called "spatial dither", introduces me ail. 
possibility of modulating or smearing "smart" jammer 
signals, thereby limiting their effectiveness. 
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APPENDIX A 



ANALYTICAL SOLUTION OF THE A.N.C. 

Consider the two-weight noise canceller of Figure 2.3, 
and let us assume that the samoled reference mo its are: 



x. T = C cos ( j ' T 
1 J o J 



X 2J = C sln( V T ‘ 51 



Figure A. 1 is a flow diagram showing signal propagac. . _ n 
s uch an adaptive noise canceller. The first step to rue- 
analysis is to ootain the isolated impulse response frit t. 
error, e , point C, to the filter output, point G, with rh-_- 
feedback loop from point G to point B broken. 

Let an impulse of amplitude, a, to be applied at point 
at discrete time, j = k; that is, 



e, = a« (j - k) 



( A . 3 ) 



The 5 ( j - k) is a Kronecker delta function, defined as 



« ( j - k) = 



for j = k 



other wi se . 



A. 4 
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Ficjure A.l. Siqnal Propaqat ion ihrouqh 



The impulse causes a response at point: D 



a C cos(u Q kT + t> ) for j = k 



0 otherwise. 



mien is the input scaled in amplitude oy the insta it-vie i -.-r 
/alue of at 3 = k. Tne signal flow path from po . n t E t 

point E is tna t of a digital integrator witn z •: .« n s f a •: 
function, 2u''(z-l) , and impulse response, 2-u'j-l. , wnere 
u ( 3 ) is the discrete unit step function 



u ( ] ) 




for 3 < 0 



otherwise . 



Convolving 2 uu(]-l) with 



Vu 



ilds a response at poi:r 



of 



W 1 J = 2 P a C cos(oj Q kT + 1)) (A, 7 ) 

where j - k + 1 . When the scaled and delayed step function 
is multiplied by the response at point F is obtained as 

y^j = 2 u a cos(u Q jT + 4> ) c o s ( uo^kT + «t> ) (A. 3 ) 
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where j > k + 1. The corresponding response at point J, 
obtained in a similar manner, is 

y 2J = 2 u m C 2 sin Q j T + $)sin(u Q kT + $ ) (A. 9; 

where j > k + 1. Combining the preceeding two equations 
yields one response at the filter output, point G: 

'/ 2 j = 2 ^ 1 C 2 cos ( Jq T {]—’<) ) 

= 2 u a c 2 u ( j-x-1 ) cos ( w 0 T ( 3 -k) ) (A.. 

Note that the aoove is a function only of (j-k) and is t.nus 
a time invariant impulse response proportional to me u:;u. 
impulse . 

A linear transfer function for the noise canceller may 
now be derived in the following manner. If the time, <, is 
set equal to zero, the unit impulse response of the linear 
time- invar i ant signal-flow path from point C to point G is 

y a = 2 u C 2 u( j-1) cos( u 0 jT) (A. 11) 

and the transfer function of this path is 

Z ( Z - COS 0) T) 

"2 2 - 1 

z - 2z cos a) T + 1 



G (z) =2 u C 2 



56 



(A. 12) 



2 

2yC (z cos u q t - 1) 

__ 

z - 2z cos jj q T + 1 



This function can be expressed in terms of a radian sampling 
frequency, 2 = 2ir/T, as 

2yC 2 [z cos(2tx a - ") - 1] 

G(2) = ' ii 

2 - 2 2 COS ( 2 T J 7 ) +1 

o 



If the feedback loop from point G to point 3 is now closed, 
the transfer function, H(z) , from the primary input, paint 
A, to the noise canceller output, point C, can oe Detained 
from the feedback formula 



H ( z) = 



2z cos ( 2 iuj 2 
o 



_ i 



+ 1 



z 2 - 2(1 - y C 2 ) z cos(2iru « L ) + 1 - 2yC 2 



(A. 14 



This equation shows that the noise canceller with a cosine 
reference input has the properties of a notch filter at the 
reference frequency along the signal flow path from primary 
input to output . 
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APPENDIX B 



TAPPED DELAY LINE FILTER 



The tapped delay line [Ref. 5] is the preprocessor .tost 

often used in adaptive filtering applications. Most 

mathematical analyses of adaptive filters are cased on 

tapped delay line model. A tapped delay line can be 

implemented either in analog or digital form. The sc uem 

of an adjustaole filter based on an analog tapped delay L : 

is shown in Figure 3.1. Such a delay line has an analog 

input and the output at each tap is a delayed version of 

t In 

this input. Thus the n 1 component of X(t) is the output at 
t n 

the n * tap and is given by 



x n ( t ) = x ( t - (n - I)*) 



(3.1) 



where & is the delay between taps and where 



X (t) 



x x (t) 
x 2 ( t) 





( B . 2 ) 
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Figure 3.?.. Tapped Delay Line Filter 
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In digital implementation, the tapped d.elay line is 
implemented as a shift register either in hardware, Figur 



B . 2 , or by simulation inside the computer, Figure 3. 
new digital sample of the filter input is taken at ev 
adaption cycle and shifted into the leftmost position 
shift register. Old samples are shifted right by one 
position and the oldest sample is shifted out altoget 
Each position of the shift register is multiplied by 
corresponding weight. Thus the digital output of to, 
is the weighted sum of the last N samples of tne moj 



e r y 
o £ 



60 



CO 



Arid 1 og 
Ci rcui ts 




o 

on 1 
oo 
O LU 

O o I 
— J O i 

< q: 

^ Cl. 

< UJ 
CH 

Cl. 



0> 

00 

5 

Cl 

00 

CD 

s- 



(/I 

CU 



3 

Cl 



cn 

O 

*3 

c 

<T3 



61 



Figure B.2. Digital Filter Type Adaptive 
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des i red res|)onse 



APPENDIX C 

COMPUTER SOFTWARE UTILIZED 



All plots through this work were obtained us in 
and CSM9 language programs running on the IB.'-'. 30 53 rain 
fra me ccaputer. 

The following CSX? language program „as use: (::: 

various different data) to prod ace the plots 2.5, c, 7, 
and 9. As inputs to the prograi are entered tne ccef f _c:-. r: ' 
of the transfer function H (z) = I { Zj /•: ( z) 5.r. t..e le.: 

Z-iomain and as output is the plot of the nagr.it uc>= versus 
freguencv of the transfer function H ( z) . 

//MANIKA1 JCB (25o7, 1236) , • K AS IS AS’ ,CLASS=E 
//* FC 5 M A Z 9 5,0311 AM E= P IC 1 . 5 Y £ VECT E , 3 ES 1=102 A L 



//X . CCMPFINT DD COilM Y 
// X. SYSPFINI DD DUMMY 
//X . PICIFAFM DD * 

EPIC! SCAI E=. 5 SEND 
//X . SYSIN DD * 

*t#*iii*4**ii**<!**^ a********************* i + 

*THI S PRCGFAM COMPOTES THE FREQUENCY RESPONSE OF UP IC A * 

* TENTH ORDER FILTER WHOSE TRANSFER FUNCTION * 

* CAN BE EXPRESSED AS: * 

* * 



// EXEC C S M P XV 
//X.SYSLI3 CD 



// 

// 

// 

// 



DD 



ED 



ED 
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H ( Z) = Y ( Z) / X (2 )=H(Z)/D (I) 



* 

* 

* 

* 

a 

* 

* 

* 

* 

* 

* 

* 

* 



l i f. w 



4c 



-1 -2 -3 -4 -10 * 

N(Z)=A0 + A1*Z +A2*Z + A3*Z +A4*Z + + A10*Z , A N C - 

-1 -2 -3 -4 -10 * 

D (Z ) = 30 + 3 1*Z +32*2 +33* Z +34*Z + ... +310*Z . * 



MAGNITUDE AMD PHASE ARE GIVEN F03 THE INTERVAL 
0 < N < ? I ONLY SI MCE THAT INFORMATION 15 
REPEATED ON THE INTERVAL FI < * <2*?I. 

PHASE IS CONSTRAINED TO PLUS OR MIN'JS loO DEGREES . 
THE 0 r D I NAT E IS IN UNITS Or <> / PI . 






INITIAL 

RENAME 1 1 M Z = AO V R FI 

FAR AM F 1 = 3 . 1415 5 

PARA M M = 0 . 0 0 1 

PAR AN SJ 1 = 0. 25 

C 1 = — 2 .*CCS ( 2. * ? I *W 1 ) 

C2=-2.*(1.-M)*C0S(2.*FI*K1) 

C3= 1.-2. *M 

A 1 =C 1 

E 1 = C 2 

£2=C 3 



* 


ENTER 


T HE 1 A ' 


COEFFICIENTS 


Or 


PARAM 


A 0 = 1 


. 0 , A 2 = 


1. „ A3= 0. , A4= 0., A5 = 


0 .* 


PAR AN 


A 6 = 


0. ,A7 = 


0 . , A 8= 0.,A9= 0 . , A 1 0 = 


0 . * 


##**#****** *** 


************* ************************ 


****** 


* 


ENTER 


T HE * E * 


COEFFICIENTS . 


4c 


PARAM 


E C = 1 . C 


,33 = 


0. , E4= 0. , 35= 0 . 


4c 


PARAM 


E6 = 


0 . ,37 = 


0. , B3= 0 . , 39 = 0 . , 3 1 0 = 


0 .* 
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DYNAMIC 

NOS C 5 1 

N=WO VRPI*PI 

DZGR Z 1 = 160. *SCVE?I 

•MAG 1 = AO + A 1 * COS (K) + A2*COS (!’* 2 . ) +A 3 *C 05 (W * 3 . ) . . . 

+A4*C0S (W*4.) +A5*CGS (fi*5.) +A6*C05 (4*6. ) 

+ A7* COS (W*7. ) +A 8 *CCS (W* 8 .) 

+ A9* COS (W*9.) + A 1 0*CO S (Tv*10.) 

MAG2=A1 *SIN (W) + A 2*SI M (W*2. ) +A3*SIN (W*3 . ) ... 

+A4*SIN (W*4.) + A5*SIN(N*5. ) + At*SIN (4*6.) .. 

+A7*SIN(N*7.)+AS*SIN(W*8.) ... 

+ A 9*SI N (.» *9 .) + A 1C* SIN ( W * 1 0 . ) 

MAG3 = rO + Z 1*C0S (W ) + 32 *COS (W* 2. ) +33*CCS (« *3 . ) . . . 

+ £4* CCS (IT *4 . ) +£5*C0S (K*5. ) +3 6 * COS (w*c . ) 

+ £ 7 * CCS ( w *7 . ) + £8 *CCS (W* 8 .) 

+ E9 * CCS (N*9. ) +E 1 0*CO S (S*10.) 

MAG4=E1*SIN (ft) +3 2* SIN (ft* 2. ) + 3 3 * S I N ( ft * 3 . ) ... 

+ 3 4 *SI N (ft * 4 .) + 55*SIN(ft *5. ) +36* SIN ( 4 * 6 .) .. 

+ 3 7 * S I N (7*7 .) + 3 £ * S I N ( ft * 8 . ) ... 

+ 3 9*SIN (ft *9 .) + 3 10* SIM ( ft*1 0. ) 

MAG=SQET (M AG 1* * 2 +M AG 2**2 ) /S QH7 ( M A G 3 * * 2 + M A G 4 * * 2) 
IE (MAGI. EQ. 0.0.AND.MAG2.EQ. 0.0) A ft G 1 = 0. 

I? (MAGI - NE.O.O.OE- M A G2 - N E. 0 .0) . . . 

ANG1=ATAN2 ( M AG2 , M AG 1 ) *180. /FI 
IF(MAG3.ZQ.0.0.AND.MAG4.EQ. 0.0) ANG 2=0 . 

IF (MAG3.NE.0.0. OH.MAG4. NE. 0 .0) . . . 

ANG2=ATAN2 (M AG4, MAC-3) *18 0./ FI 
FHASE = ANG1- ANG2 

IF (PHASE. GT. 180.) P HA SE=P H A SS-3 6 0. 

IF (PHASE. IE. -180.) PEASE=PH ASZ +3 60. 

TIMZF FINT IM=1 . 0 

PRIM ft, CEGF.EE , M AG , PEASE 

CUT? U I HC 7 RPI , MAG 

IAS E I DIGITAL FREQUENCY RESPONSE - MAGNITUDE 
IA3EI CF THE NOTCH FILTER ( KOVR?I=W/PI, M = 0 . 00 1) 
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EAGE XYEICI 
OUTPUT KCVEPI, PHASE 

IA3EI FREQUENCY RESEC NSE - PHASE IN DEGREES 
IA3EI CF I HE NOTCH FIITEE ( KOVEP I=w/?I, :i = 0 . 00 1) 

PAGE XYEICT 
END 
STOP 
END JOE 
/* 

/* 

/* 

/* 

The following FORTRAN program was used to proiace t 
tea a p at tern of figure 3.5 . 

C 

C THIS ERCGRAM DOES PCLAR PLOTTING OF A GIVEN POLA? 

C FUNCTION (VTHE(I)) .VS. THE ANGLE ( ' 7 ( I ) ) 

INTEGER I , J 

REAL K (73) , H A G 1 (73) , M AG 2 (73) , 

1 HAG 3 (73) , HAG 4 (73 ) , VTHS (7 3) , 

2 A0/A1/A2 / 30,31,32/01/02,03, 

3 MAX V , ?I , H , li 1 
C 

PI = 3. 14159 
M = . 0 1 
W 1 = . 1 25 

C1=-2.*COS ( 2 - * P I * K 1 ) 

C2 = -2 .* (1.-M) *CCS (2. *EI*M1) 

C3=1. -2. *M 
A0 = 1 . 

A 1=C1 
A 2= 1 . 

E0 = 0. 

B 1 =C2 
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32=C3 

EC 61 = 1, 7 3 

r /i (I ) = (FICAT (I) -1. 0) *5.0*0. 0 17453 

6 CONTINUE 
£071=1, 73 

NAG 1 (I) =A C + A 1 *C0 S ( r « (I) ) + A 2 * C C 3 (X (I) *2. ) 
MAG 2 (I) = A 1*SIN (W (I) ) + A 2 * S I N (X (I) *2 .) 

NAG 3 (I) =B 0+3 1 *C0 S ( t a (I) )+32*C03 (W (I) *2. ) 
MAG 4 (I) =3 1 * SIN ( 57 (I) ) +52*SIN ( N (I) *2 .) 

7 CONTINUE 

EC 10 I = 1, 73 

VTHE (I) = S £RT (MAG 1 (I) **2+MA32 (I) **2)/ 
x SQPT (MAG3 (I) **2 + HAG - (I) **2) 

< (I) = W (I) * 57 . 256 

10 CONTINUE 

C *** NCFNAIIZE AMD E CRM THE OUTPUT ARRAY VTHE (I) 

C 

MAX V=VTH E ( 1 ) 

EC 12 I = 2, 73 

IF ( (VTHE (I) ) . LS . (M AXV) ) GO TC 12 
MAXV= VTHE (I) 

12 CONTINUE 

DC 13 I = 1, 73 

VTHE (I) = VTHE (I ) /M AXV 

13 CONTINUE 
C SET UP TEK618 

CALI TEK618 
C CAII COMPRS 

CALI PAGE (8. 5,11.) 

CAII BLO WOP ( .6) 

CAII NOBEDR 

CAII AREA2D (7. 5,7. 5) 

CAII HEAEIN (’EROADSIDE ARRAY PATTERN' ,-23 
C SET UP POIAR GRID 

CAII POLAR (0.01745, . 26 6,3.7 5, 3.75) 
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CALI MARKER ( 16) 

CALI DOT 

CALI CURVE (W ,VTHE,7 3, 1 ) 
CALI RESET ( ’ EOT') 

CALL 3 BID (2 , 4) 

CALL SNDPL {0 ) 

SLOP 

END 
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